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Abstract. We study the geometry of oriented right-angled hexagons in H 4 , 
the hyperbolic 4-space, via Clifford numbers or quaternions. We show how to 
augment alternate sides of such a hexagon so that for the non-augmented sides, 
we can define quaternion half side-lengths whose angular parts are obtained 
from half the Euler angles associated to a certain orientation-preserving isom- 
etry of the Euclidean 3-space. This generalizes the complex half side-lengths 
of oriented right-angled hexagons in H 3 . We also define appropriate complex 
half side-lengths for the augmented sides of the hexagon. We further explain 
how to geometrically read off the quaternion half side-lengths for a given ori- 
ented, augmented, right-angled hexagon in H 4 . Our main result is a set of 
generalized Delambre-Gauss formulas for oriented, augmented, right-angled 
hexagons in H 4 , involving the quaternion half side-lengths and the complex 
half side-lengths. We also show in the appendix how the same method gives 
Delambre-Gauss formulas for oriented right-angled hexagons in H 3 , from which 
the well-known sine and cosine laws can be deduced. These formulas generalize 
the classical Delambre-Gauss formulas for spherical/hyperbolic triangles. 



1. Introduction 

Teichmiiller theory and the study of right-angled hexagons in the hyperbolic 
plane H 2 are closely related. In particular, in studying the Weil-Petersson metric 
and the Fenchel-Nielsen coordinates for the Teichmiiller space, a hyperbolic sur- 
face is decomposed into pairs of pants with geodesic boundary, which are further 
decomposed into planar right-angled hexagons. The side-lengths of the resulting 
right-angled hexagons then form part of the parameters of the Teichmiiller space. 
This can be extended to the study of quasi-Fuchsian space using complex Fenchel- 
Nielsen coordinates in an analogous way, using right-angled hexagons in the hyper- 
bolic space H 3 and complex side-lengths instead; see, for example, Tan [2(>], Scries 
[25], and Goldman [15]. Parker and Platis [23] have also made analogous studies 
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on representations of surface groups into the isometry group of complex hyperbolic 
plane. 

Fundamental to this point of view is that, generically, there is a direct corre- 
spondence between equivalence classes of irreducible representations p of a rank 
two free group F2 — (X,Y) into Isom + (H 2 ) (respectively, Isom + (H 3 )) with points 
in the moduli space of right-angled hexagons in H 2 (respectively, H 3 ): the axes of 
p(X), p(Y) and p(Y~ 1 X~ 1 ) form the alternate sides of a right-angled hexagon, 
whose (complex) lengths are half the (complex) translation lengths of p(X), p(Y) 
and p(Y~ 1 X~ 1 ), respectively. 

Things are more complicated for representations of F 2 into Isom + (H 4 ), where 
H 4 denotes the hyperbolic 4-space. Although one can still define a right-angled 
hexagon from a generic representation p as above, the equivalence class of p is not 
determined by the right-angled hexagon, which can be seen by a simple dimension 
count. So there is a non-trivial fibre for the map from the character variety 

Hom(F 2 , Isom+ (H 4 ))/ /Isom+ (H 4 ) 

onto the space of equivalence classes of right-angled hexagons in H 4 . Thus, the 
study of this character variety is significantly more difficult and we hope to pursue 
this in the future. As a preparation for this, this paper is devoted to studying 
the geometry of right-angled hexagons, in particular, to obtaining trigonometric 
relations such as the law of cosines or its equivalent for the geometric configuration 
determined by p(X), p(Y) and p{Y- 1 X~ 1 ), 

Recall that a convex right-angled hexagon in the hyperbolic plane H 2 is a cycli- 
cally ordered sextuple (Li, ■ ■ • , Lq) of complete geodesies in H 2 such that for each 
n = 1, ••• ,6, the two geodesies L n and L n+ \ (with indices modulo 6) intersect 
perpendicularly (and these are the only intersections among the six geodesies). 

It is well known (see, for example, Beardon's book [ ], pp. 160-161) that, for a 
convex right-angled hexagon (L\, • • • , Lq) in H 2 with side-length l n > for the side 
L n , n = 1, • • • ,6, the Law of Cosines holds: for each n modulo 6, 

cosh/ n = — cosh l n +2 cosh l n +4 + sinhZ n+ 2 sinhZ n+ 4 cosh/„ + 3, (1) 

as well as the Law of Sines: 

sinhZi sinh/3 sinh^s . 
sinh Z4 sinh Iq sinh I2 

An oriented right-angled hexagon in the hyperbolic 3-space H 3 is a cyclically 
ordered sextuple (L 1: ■ ■ ■ ,L 6 ) of oriented, complete geodesies in H 3 such that for 
each n modulo 6, the two geodesies L n and L n+ \ intersect perpendicularly. 

Each ordered triple (Z n _i, L n , L n +i) defines a complex length <r n for the side 
L n , defined up to multiples of 2ni, that is, a n G C/27rzZ. This is the signed complex 
translation length of the Mobius transformation with oriented axis L n which sends 
L n -i to L n +i (see §7.3 for the precise definition). 

By Fenchel's work (see [14, Chapter VI]), the complex side-lengths {a n }n=i °f 
an oriented right-angled hexagon in H 3 satisfy several trigonometric identities, the 
most important of which are the Law of Cosines: for each n modulo 6, 

cosh er n = cosh er n +2 cosh <7 n +4 + sinh <J n -\-% sinh cr n +4 cosh cr n +3, (3) 

and the Law of Sines: 

sinhci sinh 03 sinh (75 . . 

sinner sinh (76 sinh<72 
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(Note that the above law of cosines and that of sines for oriented right-angled 
hexagons in H 3 were known to Schilling as early as in 1891 [24], but the correct 
treatment of signs and a convincing proof seem to be first given by Fenchel in [14].) 

One nice aspect of the law of cosines and the law of sines above for oriented 
right-angled hexagons in H 3 is that they specialize to various well-known identities 
for relatively simple polygons, convex or self-intersecting, in H 2 , such as triangles, 
quadrilaterals with two right angles, pentagons with four right angles, and right- 
angled hexagons; see [14, §VI.3] for a complete list of formulas in each case. In fact, 
this can be seen by observing that each of these planar polygons can be regarded as 
a subset of the sides of a right-angled hexagon in H 3 ; for example, the three sides of 
a triangle may be regarded as three alternate sides of a right-angled hexagon in H 3 . 
In particular, in the case of a convex, planar right-angled hexagon, we may regard 
it as an oriented right-angled hexagon in H 3 by orienting its side-lines consistently. 
The complex side- lengths are then l n + ni, n = 1, • • • ,6 and formulas (3) and (4) 
above simplify to formulas (1) and (2). 

On the other hand, recall that, for a spherical triangle in the unit sphere, there 
are the Delambre-Gauss formulas (142)-(145) (often called Delambre's analogies 
or Gauss formulas separately) for the half side-lengths and half interior angles of 
the given triangle, from which the law of cosines and the law of sines follow easily. 
Similar formulas hold for hyperbolic triangles in the hyperbolic plane. 

In §7, we obtain Delambre-Gauss formulas for oriented right-angled hexagons in 
H 3 (Theorem 7.9), from which one may easily obtain (3) and (4), namely the law 
of cosines and the law of sines, for oriented right-angled hexagons in H 3 . 

The main purpose of this paper is to obtain analogous formulas for oriented 
right-angled hexagons in H 4 . For convenience, we shall use the simple terms lines 
and planes to mean respectively complete geodesic lines and totally geodesic 2- 
planes in H 4 . To start with, we observe that for a given oriented right-angled 
hexagon (Li, • • • , L§) in H , one could similarly define a "complex length" for each 
side L n , which is invariant under orientation-preserving isometries of H 4 . The six 
complex lengths so defined, however, do not determine the right-angled hexagon up 
to isometries of H , as the information is incomplete. We overcome this problem by 
considering oriented, augmented, right-angled hexagons as follows. Starting with 
an oriented right-angled hexagon (L\, . . . , Lq) in H 4 , we choose three alternate side- 
lines, say L/2, £4, Z/6) and augment them as oriented line-plane flags F n = (L n ,II n ), 
n = 2,4, 6, where F n consists of L n together with a oriented plane n n in H 4 such 
that L n is contained in II„, and, crucially, the lines and L n+ \ are both 

perpendicular to II„; see Figure 1 for an illustration. 

It is always possible to define an oriented, augmented right-angled hexagon from 
a given oriented right-angled hexagon in H 4 . In fact, generically, the second condi- 
tion above ensures that II n is well-defined up to orientation, unless L n -i,L n and 
L n+ i are coplanar, in which case there is an S 1 worth (or more) of possible II„'s. 
We may of course alternatively augment the three odd-indexed side-lines. However, 
we do not consider a full augmentation of all the side-lines. 

Consider an oriented, augmented right-angled hexagon (L%, F2, L3, F4, L5, Fq). 
Let A2 be the Clifford algebra Clo,2, that is, the associative algebra over the reals 
generated by ei,e2, subject to the relations eie2 + e2ei = and e\ = e\ = — 1. As 
ungraded algebras, A2 is isomorphic to the algebra of quaternions. 
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Figure 1. An oriented augmented right-angled hexagon in H 4 

For each side-line L n , n = 1,3,5, the configuration (F n —i,L n ,F n +i) allows us 
to define two quaternion half side-lengths 5 n € A 2 modulo period, where the two 
values of 5 n modulo period differ by iru for some u G A2 with u 2 = — 1. They 
can be described by a real translation part and three angles. The real part is 
uniquely defined to be half the signed translation length of moving the intersection 
point of L n with L„_i to that of L n with measured along the oriented axis 

L n . The three angles are determined by the orientation-preserving isometry of the 
Euclidean 3-space R 3 = R + Rei + Re 2 induced by parallel translation along L n 
which sends the unit tangent vector on the unit sphere S 2 in R 3 determined by the 
oriented flag F n -i to that determined by F n+ i; explicitly, they are the half Eulcr 
angles associated to the above isometry of R 3 . On the other hand, for each flag F n , 
n = 2,4,6, the configuration (£ n _i, F n , £ n +i) allows us to define two e2-complex 
half side-lengths modulo 27re2, denoted by S n , that is, 5 n € (R + e2R)/27re2Z. Note 
that the two values of S n modulo 27re2 differ by 7re2- 

The main result of this paper are the generalized Dclambre-Gauss formulas, (5)- 
(8) below, for oriented, augmented, right-angled hexagons in H 4 . We expect that 
these formulas will prove useful in the study of hyperbolic 4-manifolds as well as 
the representation varieties of surface groups into Isom + (H 4 ). 

We remark that the asterisk notation ()* appearing in formulas (5)-(8) denotes 
the reverse involution of the graded algebra A2 (see §2.2): 

(.t + x x e\ + x 2 e 2 + xi 2 eie 2 )* := a; + x x e\ + x 2 e 2 - x 12 e x e 2 , 

with real coefficients Xo, x\, x 2 , x\ 2 , while the definitions of the hyperbolic functions 
cosh and sinh with an A2-variable will be discussed in §3.3: in short, we define 

exp(x) + exp(— x*) . exp(x) — exp(— x*) 

coshx := , smlix := . 

2 2 

Theorem 1.1 (Generalized Delambre-Gauss formulas for hexagons in H 4 ). For an 

oriented, augmented, right-angled hexagon (Li, F2, L3, F4, L5, Fq) in H with any 
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choice of {ei, e-i\ -quaternion half side-lengths 61,63,65 and e2-complex half side- 
lengths 62,64,66, the following formulas hold: 

sinh 61 cosh 62 sinh 83 + cosh 6± cosh 5% cosh 63 

= e (sinh 84 cosh £5 sinh 5§ + cosh 64 cosh 65 cosh 6$)*; (5) 
sinh <5i cosh 62 cosh + cosh <5i cosh 82 sinh £3 

= e (sinh (5 4 sinh 85 sinh <5 6 — cosh £4 sinh 85 cosh 86)*; (6) 
sinh £1 sinh 82 sinh £3 — cosh 81 sinh 82 cosh ^3 

= e (sinh (5 4 cosh 85 cosh 5 6 + cosh 84 cosh 85 sinh <5 6 )*; (7) 
sinh 8 1 sinh 62 cosh ^3 — cosh 8 1 sinh 62 sinh ^3 

= e (sinh (5 4 sinh 85 cosh 56 — cosh 84 sinh <5s sinh 6q)* , (8) 

wit/i £ = 1 or — 1, depending on the choices of the six half side-lengths {<5 r i}n = i- 

Remark 1.2. We emphasize that, for a fixed choice of the six half side-lengths, it 
is the same £ that occurs in the identities. For easy memorization, the complicated 
formulas (5)-(8) above can be conveniently abbreviated as follows: 

(scs + ccc)i23 = £(scs + ccc)4 56 ; 
(scc + ccs)i23 = ^(sss - csc)4 56 ; 
(sss-csc)i23 = £(scc + ccs)4 56 ; 

(SSC-CSS)123 = £(SSC - CSS)4 56 . 

Remark 1.3. These formulas rewritten in another form, (138)-(141), which makes 
use of the operations © and G in A2 to be defined later will be given in §6. 

To fully understand the meaning of these identities, and the terms involved, it 
is useful to understand the difficulties involved in attempting this generalization. 
The main difficulty is related to the fact that the point-wise stabilizer of a com- 
plete geodesic in Isom + (H 4 ) is isomorphic to SO(3), which is non-commutative, as 
opposed to the case of H 3 , where the point- wise stabilizer of a complete geodesic 
in Isom + (H 3 ) is isomorphic to SO(2) and is commutative. To a certain extent, 
this non-commutativity is reflected by the non-commutativity of A2, so that the 
representation of elements of Isom + (H 4 ) by Vahlen matrices (2x2 matrices with 
entries in A2 satisfying certain conditions, following Vahlen [27] and more recently 
Ahlfors and his collaborators and several others ([1]-[6],[21], [28] ,[29] ,[30])), is par- 
ticularly useful and appropriate. This is the approach adopted here. It also has the 
advantage of shedding some light on the geometry of A2, in relation to H 4 and 9H 4 ; 
indeed some of the results in §3 and §4 are of independent interest, in particular, 
the definitions of the functions exp, sinh, cosh and log for elements of A2 and the 
Euler decomposition of units of A2 (Proposition 4.1 and (91)). Nonetheless, this 
lack of commutativity also means that appropriate generalizations of the hyperbolic 
functions sinh and cosh to Clifford numbers (or quaternions) need to take careful 
account of the non-commutativity. We do this in §3.3 and it is with respect to 
these definitions that the identities (5)-(8) in Theorem 1.1 should be interpreted. 
The non-commutativity also means that we are unable to obtain analogues of the 
cosine and sine rules from the Delambre- Gauss formulas in H 4 . 

The rest of this paper is organized as follows. In §2 we introduce the Clifford alge- 
bra A n = Clo,n and briefly review the theory of representing Mobius transformations 
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of general dimension by 2 x 2 Vahlen matrices. In §3 we focus on the Clifford algebra 
A2, and carefully define the exponential function exp, the hyperbolic trigonometric 
functions cosh and sinh, the logarithmic function log of an A2-variable, as well as 
the operations © and © in A2. We also give the Vahlen matrix representations 
of Mobius transformations of <9H 4 which fix —1 and 1 in terms of the hyperbolic 
trigonometric functions defined. In §4 we study the geometry of SO (3) via the 
Clifford algebra A2; in particular, we study the Euler decomposition of units in A2 
and explore the algebraic and geometric meaning of the associated Euler angles. 
In §5 we define {ei, e2j-quaternion half side-lengths and the e%- and e2-complex 
ones for configurations (F n -\, L n , F n +i) and [L n -i, F n , L n +i) of oriented lines and 
flags in H 4 . In §6 we prove our main theorem (Theorem 6.5, rephrasing Theorem 
1.1), the generalized Delambre- Gauss formulas involving the {e±, e2j-quaternion 
and e2-complex half side-lengths for oriented, augmented, right-angled hexagons in 
H 4 . Finally, in the appendix (§7) we establish the Delambre-Gauss formulas for 
oriented right-angled hexagons in H 3 , whose proof uses the same method used in 
the proof of Theorem 6.5 but does not require the material in the earlier sections. 
The reader who is familiar with the geometry of H 3 and would like to understand 
the geometric idea of the proof of the main theorem should first read the appendix. 

Acknowledgements. The authors would like to thank Bill Goldman, Roger Howe, 
Sadayoshi Kojima, Frangois Labourie, John Parker, Makoto Sakuma, Caroline Se- 
ries, Masaaki Wada, Hongyu Wang, Yanlin Yu, and Qing Zhou for helpful conver- 
sations. 

2. The classical Clifford algebras A„ 

2.1. The classical Clifford algebras A n . Following Ahlfors, we use A„ to denote 
the classical Clifford algebra Clo, n , that is, the associative algebra over the reals 
generated by the elements ej., e%, ■ ■ ■ , e„ subject to the relations 

2 _ 2 _ _ 2 _ 1 

e l — e 2 — ' ' ' — e n — L > 

CjCj + CiCj = for i 7^ j. 

Thus A„ is a real algebra of dimension 2™ and is a subalgebra of A n +i. 
An element of A„ can be written uniquely in the form 

^ <^0 © ^ dii^ix © ^ ^ j 2 i 3 6i x C^ 2 ej 3 © 

l<zi<n l<ii<Z2<n l<ii<Z2<*3<n 

+ + ai2-n eie2 ■ • • e„ (9) 

with real coefficients a , a, l: a ili2 , a lli2 i 3 , • • • , a ili2 ... ln _ 1 , a 12 ... n - 
Let us write for the degree p part of a, that is, = ao and 

a (p) = a ili2 ... lp e ll e l2 ---e ip , p = l,- ■■,<!%. (10) 

l<ii <i2 <- -<i p <n 

Then 

a = + + + + ■■■+ fl W. (11) 

We call the elements in A„ of the form the (degree) p-vectors in A„; in 
particular, the 0-vectors are the reals. Correspondingly, we have a decomposition 
of A„ as the direct sum of its p- vector subspaces (with = K) : 

A„ = © A^ © • • • © A^ . (12) 
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Notice that Ao = R, Ai = C, and A2 = H, the quaternions. 

An element of A n is called even if it is the linear combination of even-degree 
vectors of A n ; similarly for odd elements. The even elements of A„ form a subalgebra 
A+ of A n , while the odd elements only form a subspace A~ of A n . 

We shall conveniently regard A n as a Euclidean space with norm 

M = (El fl 'M 2 -J 2 ) 1/2 (13) 

(where the sum runs over all its coefficients) and the inner product 

(a, b) = a ili2 ... lp b lll2 ... ip . (14) 

2.2. Three involutions of A„. There are three involutions of A„. The main (or 
prime) involution ()' : A„ — > A n , the reverse (or star) involution ()* : A„ — > A n and 
the conjugate (or bar) involution () : A n — > A„ are defined respectively by 

a' = a ^ -a^ -a^ +--- + (-l) n a^; (15) 

a* = a^+a^-a^-a^ + --- + (-l) n ^- 1 ^ 2 a^; (16) 

a = a (°) -a« -a^ + a® + ■ ■ ■ + (-l) n(n+1 ^ 2 a (n K (17) 

In particular, for the basis elements, we have 

(ei ± e l2 ■ ■ ■ e ip Y = (-e ll )(-e !2 ) • • • (-e ip ) = (-l) p e u e i2 • • • e lp ; (18) 

(Cj'i 6j 2 ' ' ' e i p ) = e ip6i p _ 1 £i 1 ( — 1) P ^ P ^ Bjj &i 2 ■ ■ ■ &i p ; (19) 
Cil 6i 2 ' ' ' e,; p = ((Ci! 6j 2 ' ' ' e i P ) ) ( 1) P ^ P ^ 6ix 6j 2 ' ' ' e *p ■ (^0) 

Remark 2.1. Note that in [20] the main involution ()' is denoted by ()" and the 
reverse involution ()* is denoted by ()~. 

It is easily checked that the composite (in any order) of any two of the three 
involutions above is the remaining one, that is, 

(a')* = (a*)' = a, W=(a)*=a', (a)'=W=a*. (21) 

Note also that the main involution is an algebra isomorphism and the other two 
involutions are anti-isomorphisms; in other words, they are all isomorphisms for 
addition, and, for multiplication, 

(ab)' = a'b', (ab)* =b*a*, ab = ba. (22) 

For the rules of commutativity, we note the following special cases: 

aei = aa if a does not contain a; (23) 

aei = —&io! if all terms of a contain e,. (24) 

2.3. The group A* of invertible elements in A„. As usual, an element a G A„ 
is said to be invertible if there exists an element b € A n such that 

ab = ba = 1. (25) 

By E. Cartan [10] , A n is isomorphic to a subalgebra of the matrix algebra Mat(m, M) 
for some integer m; hence the single equality ab = 1 implies ba = 1. Such an element 
b € A„, if it exists, is unique and is denoted by a -1 , as usual. The set of all invertible 
elements of A n is a multiplicative group which we denote by A^ . 
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As simple examples, it is easy to verify that 

A* = R x = R\{0}, A? = Ai\{0}, A 2 X = A 2 \{0}, (26) 
A 3 X = A 3 \((l - eie 2 e 3 )R U (1 + eie 2 e 3 )E). (27) 

2.4. The spaces of 1-vectors and para- vectors of A„. We are interested in 
the 1-vector subspace A^ and the para-vector subspace A^ 0,1 ^ := A^ A^ which 
are of (real) dimensions n and n + 1 , respectively. 

It is clear that all the three involutions leave the 1-vector subspace A„^ and the 
para-vector subspace A^ ' 1 ^ invariant. 

Lemma 2.2. For x G An , we have x* = x, x = x' , and xx = xx — \x\ 2 , while 
for x G An \ we have x = —x and x 2 = — \x\ 2 . 

Hence every non-zero para-vector x G An ' 1 ^ is invertible, with x" 1 — i/|a;| 2 . 
As for the inner product in An = M n+1 , we have 

2{x,y) =xy + yx for x, y e A^ ' 1 ). (28) 

2.5. The pure Clifford and full Clifford groups of A„. Following Chevalley 
[12], we call the multiplicative group consisting of all products of non-zero 1-vectors 
in A„ the (pure) Clifford group of A„ and denote it by rP urc ; that is 

rr = {aeA n \ a = x lX2 ---x m ^Q, X!,x 2 ,--- , x m G A^ 1 ', m > 0}, (29) 

where, as a convention, we regard the null product as 1. (Note that the idea of the 
definition goes back to R. Lipschitz [18].) Following Ahlfors-Lounesto [6], we call 
the group consisting of all products of non-zero para- vectors in A n the full Clifford 
group of A„ and denote it by T^ 1 ' 1 or simply r„; that is, 

r„ = {a G A„ | a = xix 2 ■ ■ ■ x m ^ 0, Xi,x 2 ,- ■■ , x m G A^ ' 1 ), m > 0}. (30) 
Notice that rP ure is subgroup of r n . As simple examples, we have 

r = A x , ri = Af, r 2 = a 2 x = a 2 \{o}. (31) 

We also define the pure even Clifford group r+ purc of A n by 

p+pure = r Purc nA + 

= {a G A„ | a = xix 2 ■ ■ ■ x 2m ^ ° 5 %i,X2, ■ ■■ , x 2m G A^ 1 ), m > 0}. (32) 
The two versions of Clifford groups are related by the following proposition. 

Proposition 2.3. r+£" re r„. 

Proof. Notice that A„ = A^ +1 as real algebras with an isomorphism given by 
e, <-> eie„+i, i = 1, 2, • • • ,n. The desired isomorphism follows by observing that 

(siei H h s n e n + s„ + ie„ + i)(tiei H h t n e n + t n+ ie n+ i) 

= (sieie„ + i H h s„e„e„+i - s n+ i)(*ieie n+ i H h t n e n e n+ i + t n+ i), 

with real coefficients Sj and tj, i = 1, 2, • • • , n + 1. □ 

Proposition 2.4 (Waterman [30]). (i) For a G r„, aa — aa = \a\ 2 G M x . 
(ii) For a,b G A n with a G T n or b G r„, \ab\ — |a||6|. 

The following deep theorem says that an clement of T n with real part 1 is actually 
determined by its 1- and 2-parts; for a proof, see [19] or [6]. 
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Theorem 2.5 (Lipschitz- Vahlen) . Given arbitrary u G A„ and v £ A^ 2 - 1 , £/iere 
exists exactly one element a G r„ with a*- ' = 1, a^ — u and = v. 

As a consequence of Theorem 2.5, we conclude that the full Clifford group F„ is 
a Lie group of dimension 1 + n + n(n — l)/2 = (n 2 + n + 2)/2. On the other hand, 
it can be shown that the group A* of invertible elements in A„ is a Lie group of 
dimension 2". Hence we have T n C A* and r„ ^ A* for n > 3. 

To give an explicit example of invertible elements of A„ which are not in r„, 
consider the element a = 1 + teie2 ■ • • e n G A„ with t G R\{ — 1, 1}. Then a is 
invertible since (1 + te\e2 • ■ • e„)(l — te\e2 ■ ■ ■ e„) = lii 2 el x . That a ^ F„ (with 
n > 3 and t ^ 0) follows immediately from Theorem 2.5 since 1 G T n and a^l. 

2.6. Another characterization of the full Clifford group F„. It is easy to 
show (say, by induction) that if a G T n then for all x G a1 0,1 \ ax(a') -1 G A^ ' 1 ^. 
Conversely, we have: 

Proposition 2.6 (Vahlen-Maass). // a G A n is invertible and for all 
aa;(a')" 1 G Al ' 1 ^ i/iera a G F„. □ 
Thus we obtain the following characterization of the full Clifford group: 

T„ = {a G A* I axia'y 1 G A^' 1 ) for all x G A^ 0,1 ^}. (33) 

Note that F„ is a subgroup of T„ +1 . In particular, we have, for a G F„, 

/ /\— 1 11 1 \ — 1 

ae n+ i(a ) = e n+1 a (a ) = e n+1 . 

2.7. The isomorphism SO(n + 1) = F„/1R X . An element a G F„ defines a linear 
transformation p{a) : A^ 0,1 -* — > An by 

p{a)x = axia'Y 1 for x G . (34) 

Since |p(a)x| 2 = (ax(a')" 1 )(ax(a')" 1 )' = ax(a')^ 1 aV a - 1 = |x| 2 for all x G A^, 
we see that p(a) G 0(n +1). It can be shown that p(a) G SO(n + 1) (say, first 
prove the conclusion for a = l,ei,--- ,e„, then for a G A^'^jO} by continuity, 
and finally for a G F„ by induction). This defines a homomorphism of groups 

p:T n ^SO(n + l). 

It can be shown that p is an epimorphism with kernel R x . Thus we have obtained 
the following isomorphisms of groups: 

Spin(n + 1) = T n /R+ and SO(n + 1) = F„/M x . (35) 

2.8. The one-point compactification An of A^ 0,1 -*. Following Ahlfors, we 
identify the Euclidean space R n+1 with A^ ' 1 -*, and hence the one-point compacti- 
fication K™ +1 of R" +1 with 

A! ' 1 ) := A^ 1 ' U {00}, (36) 

where, as usual, the symbol 00 := _1 operates as follows: 

oo + a = a + oo = oo for a G A„; (37) 
00 • 00 = 00; (38) 
00 • a = a • 00 = 00 for a G A* ; (39) 

X)- 1 = 0, (40) 

and the operations 00 ± 00, 00 • and • 00 are forbidden. 
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2.9. An observation of Ahlfors. Ahlfors made the following useful observation: 

Lemma 2.7 (Ahlfors [ ]). Suppose a,b £ T n . Then a~ 1 b £ An' 1 " 1 if and only if 
ab* G A.k° 4) , while aft" 1 G A-L 0,1 ' if and only if a*b £ A { n A) . □ 

2.10. The group SL(2,r rl ) of Vahlen matrices of dimension n. A Vahlen 
matrix of dimension n is a matrix ( ^ ^\ such that 

(i) a,b,c,d£ T n U {0}; 

(ii) ad* - be* = 1; 

(hi) a-^c^d^ac^.bdr 1 £ A^' 1] = a!° 4) U {oo}. 
By Ahlfors' observation (Lemma 2.7), the condition a~ 1 b £ An 0,1 ' 1 in (hi) is equiv- 
alent to ab* £ An , and similarly for the other conditions in (iii). 

On the other hand, while keeping (i) and (ii) unchanged, one can drop any two 
except the last two of the four requirements in (iii). Precisely, we have 

Proposition 2.8 (Ahlfors [2][3]). Suppose a,b,c,d £ T n U {0} and ad* ~ be* = 1. 

Then 

(a) a^,c- 1 deAl ' 1) ac-K^d- 1 £ A [ ° A) ; 

(b) a-^ac' 1 £ Ai ' 1 ) c-^^bd- 1 £ a1 0,1) ; 

(c) a-'b^d- 1 £ A^ A) c~ x d,ac~ x £ ii 0,1) . 

Example 2.9. We give a simple example to show that the converse implication of 
part (a) in Proposition 2.8 is not true, that is, 

ac^^bd- 1 £ A^ =f=> a- l b,c- x d £ A^\ (41) 

For this, let n = 2 and 

a = d=l+>§e 1 , b = c= (l - ^e 1 )e 2 £ T 2 . 

Then ad* - be* = 1, a*c = b*d = |e 2 £ A^° 4) and ab* = cd* = (| + V2ei)e 2 £ 
A 2 0,1 \ By Ahlfors' observation, ac~ 1 ,bd~ 1 £ A^, 0,1 ' and a~ 1 b,c~ 1 d (I A 2 0,1 ^ . 
Proposition 2.10. The Vahlen matrices of dimension n form a group under matrix 
multiplication, with the inverse of A — ^ ° ^\ given by 

T-) ■ <-> 

For a detailed proof of Proposition 2.10, see Waterman [30]. 

Notation 2.11. We denote by SL(2,r„) the multiplicative group of all Vahlen 
matrices of dimension n. 

Since T n c r n+ i, it follows that SL(2,r„) is a subgroup of SL(2,r n+ i). 

2.11. Mobius transformations of A^ 0,1 ' via Vahlen matrices. A Mobius trans- 
formation of A^ ' 1 ^ is defined to be a conformal automorphism of A n °' l \ or equiva- 
lently, the composition of an even number of inversions in n-spheres and reflections 
in hyperplanes in A^ ' 1 ^. In 1902, Vahlen [27] initiated the study of Mobius trans- 
formations of A^ 0,1 -* via Vahlen matrices. This study was later revived in 1949 by 
Maass [22] and re-initiated by Ahlfors in the 1980's (see [l]-[6]). 
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A Vahlen matrix A = ^ a ^ ^ of dimension n gives rise to a Mobius transfor- 
mation Ta of A^ 0,1 ' defined by 

T A {x) = {ax + b)(cx + d)- 1 £ x £ A^' 1 ); (43) 

T A (oo) = ac- 1 . (44) 

Note that Ta (oo) equals the limit of Ta{x) (where x £ A^ 0,1 ^) as x — > oo. To see 
that T A (x) £ An 0,1 ' 1 for x £ An' 1 , first note that this is the case when 



A = 



1 y\ (90 \ _ (0 -i 



11' VO g*' 1 ) \l 



where y £ An' and g £ T n . For a general A £ SL(2, T n ), this can be seen from 
the following decomposition of a general Vahlen matrix into simpler ones. 

Proposition 2.12 (Ahlfors [2] [3] [4]). A Vahlen matrix (® ^ £ SL(2,T„) can 
be decomposed into a product of simple ones as follows: 

: °T') ( c *o 1 °)(?"o)C c i d ) (45 » 

(s !)-(: °?) =»-). <«) 

Notice that —A gives rise to the same Mobius transformation as A does. Fur- 
thermore, it is easy to verify that each Mobius transformation of A^ ' 1 ^ is given by 
a Vahlen matrix and that A £ SL(2,r n ) gives rise to the identity transformation 
of A^°' 1 ' ) if and only if A = ±1. Thus a Mobius transformation is given by exactly 
two Vahlen matrices ±A and we have the isomorphism of groups 

M6b(A( l ' 1 )) S PSL(2, T„) := SL(2, T n )/{±I}. (47) 

Observe that the inclusion of groups SL(2, r n ) C SL(2, r„+i) induces the Poincare 
extension of Mobius transformations: M6b(Al°' 1 ^) C MobfA^j). 

2.12. The group Isom + (H"+ 2 ). Since the Euclidean space R n+2 is identified with 
A^+ii the upper half-space model of the hyperbolic (n + 2)-space H n+2 is given by 

n+1 

H' !+2 = (x £ A^V \ x = x + J2 x ^ x ^ x u ■ ■ ■ , Xn+i £ R, x n+1 > j, (48) 



i=l 



equipped with Riemannian metric ds 2 = {dxQ+dxf + - ■ ■ + dx 2 l+1 )/x 2 l+1 of constant 
sectional curvature —1. Its boundary at infinity, <9H™ +2 , is then identified with 
a! 0,1 -* = aI ' 1 " 1 U {oo}. It is well-known that a totally complete geodesic m-plane, 
1 < m < n + 1, is the upper half of either a Euclidean m-plane or a Euclidean m- 
sphere, both orthogonal to An 0,1 -* = M" +1 . In this model the orientation-preserving 
isometries of H" +2 are exactly the Mobius transformations of A^ 0,1 -* extended to 
H™ +2 . Thus we have the isomorphisms of groups 

Isom+(H n+2 ) S Mob^ ' 1 ') = PSL(2,T n ). (49) 
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2.13. Mobius transformations of Ai ' 1 ' fixing both and oo. 

Proposition 2.13. A Mobius transformation o/Al ' 1 ' fixes both and oo if and 
only if its Vahlen matrices ±A are of the form 

A=( a Q aer„. (50) 

By Poincare extension, A gives a Mobius transformation of A^J 1 ^ with 
7U(e n+ i) = ae n+1 a* = ade n+1 = |a| 2 e„ +1 ; 
in particular, T_4(e n +i) = e n+ i if and only if \a\ = 1. 

2.14. Mobius transformations of An 3 ' 1 ' fixing both — 1 and 1. 

Proposition 2.14. A Mobius transformation of An ^ ,/ixes 6o£/i —1 and 1 if and 

only if its Vahlen matrices ±A are of the form 

"-(; i) < 5i » 

where a, b £ T„ U {0} satisfy ab* £ A„ and aa* — 66* = 1. 

Note that ab* £ A.„ implies a*b £ An' by Proposition 2.8. However, ab* £ 
An' 1 ' does not imply a*b £ An' 1 ', as shown by Example 2.9. 

2.15. Mobius transformations of An' 1 fixing e n +i- 

Proposition 2.15. A Mobius transformation of An fixes e n+1 £ A^V^ if and 
only if its Vahlen matrices ±A are of the form 

"-(4 i) < 52 » 

where a, b £ T n U {0} safe/y a6* G a1° 4) and |a| 2 + |6| 2 = 1. 

Proposition 2.16. The subgroup o/M6b(Al 0,1 ') consisting of Mobius transforma- 
tions fixing e n +i £ A^V^ * s isomorphic to SO(n + 2). 



3. The Clifford algebra A 2 

As our aim is to study the geometry of hyperbolic 4-space, we shall be focusing on 
the case n = 2 in the rest of this paper. 



3.1. A brief account of Aa. Let us first briefly recall the basic facts that will be 
used frequently in the rest of this paper. As associative algebras over the reals, 
the Clifford algebra A 2 = Cl ,a = R + Rei + Ke 2 + Meie 2 (with e\ = e\ = -1 and 
e\e-2 + e 2 ei = 0) is isomorphic to the algebra H = E + Mi + ffij + Mk of quaternions, 
with an algebra isomorphism given by 

ei < — > i, e 2 i — yj, eie 2 < — > k. (53) 
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However, A2 is a graded algebra, equipped with three involutions; explicitly, for 
a = ao + aid + aiei + a\ie\ei G A2 with ao, ai, 02, ai2 € K, the images of a under 
the three involutions ()', ()* and () of A2 are respectively given by 

a = a - aid - a 2 e 2 + ai 2 eie 2 ; (54) 
a* = ao + aid + a2e2 — ai2de2; (55) 
a = a — aid — a 2 e 2 — ai 2 eie 2 . (56) 

An element a £ A2 is invertible if and only if a 7^ 0; in particular, a -1 = \a\~ a. 
Hence the group of invertible elements of A 2 is A 2 = A 2 \{0}. We denote the 
subgroup of A 2 consisting of all the unit elements of A2 by A™ 1 ', that is, 

A^ nit = {a G A 2 I \a\ = 1}. (57) 
The space A^ ' 1 -* of para-vectors of A2 and its one-point compactification A^ 0,1 -* are 
A^°' 1} = K + Kei + Ke 2 , A 2 0,1) = A 2 0,1) U {00}. (58) 

The full Clifford group r 2 of A2, defined as the multiplicative group consisting of 
the products of non-zero para- vectors of A2, is identical with A 2 ; that is, 

T 2 = A 2 X = A 2 \{0}. (59) 
A Vahlen matrix A G SL(2,r 2 ) is a 2 x 2 matrix ^ such that 

(i) a,6,c,de A 2 ; (ii) ad* - be* = 1; (hi) ab\ cd* G A^ 1} . (60) 
As a special rule of commutativity in A2, we notice that 

aeie2 = de2a', a G A2. (61) 

3.2. The exponential function cxp. As usual, we define the exponential function 
exp : A2 — > A2 by its Taylor expansion: 

~ x m 

expx = > — -, x G A 2 . (62) 

m! 

m— 

The convergence is guaranteed as usual since we have \x m \ = \x\ m for x G A2 and 
for integers to > 0. Note that in fact exp(a?) G A 2 since we have 

exp a; exp(— x) — exp(— x)expx = 1. (63) 

However, in general, we have, for x, y G A 2 , 

exp a; exp y ^ exp y exp a; ^ exp (a; + y) . (64) 

3.3. Hyperbolic functions cosh and sinh of an A2-variable. For reasons that 
will be clear in §3.8, we define the hyperbolic functions cosh : A 2 — > A2 and sinh : 
A 2 — >• A2, respectively, by 

coshx= e_^x+^{-x*)_ (65) 
sinhx ^expx-eM-x^ (66) 



14 



S. P. TAN, Y. L. WONG, AND Y. ZHANG 



(note the reverse or star involution in the expressions). In general, 

coshx * E ra» (67) 
00 x 2m+l 

sinh a; 7^ > — rr. (68) 

r ^ (2m + 1 v ; 

m=0 V ; 

However, the following two identities hold: 

cosh x + sinh a; = exp a;; (69) 

cosh a; — sinh a; = exp(— x*). (70) 

Since exp(x*) = (exp a;)*, cosh(x*) = (cosha;)* and sinh(a;*) = (sinha;)*, we may 
abuse notation and write them as exp x* , cosh x* and sinh x* respectively. It is 
easy to verify that 

cosh(— x) = cosha;*, (71) 
sinh(— x) = — sinh a;*, (72) 
cosha; cosha;* — sinh a; sinh a;* = 1. (73) 

3.4. The polar decomposition of non-real elements in Aa\M. It is easy to 
verify that every non-real element a G A 2 \R can be written as 

a = |a|(cos0 + usind), (74) 

where 9 G (-tt, 0) U (0, tt) and u G A^ 1,2) n A£ nit are determined by a up to sign. 
Indeed, a can be written as above in exactly two ways, the other one being 

a = |a|(cos(-0) + (-«) sin(-6»)) (75) 

with the same and u as in (74). It is useful to observe that (— &)(— u) = Ou. 
We remark that the u appearing in (74) is a square root of —1 in A2; indeed, 

A^ 2) nA™' = {u G A 2 I u (0 > = 0, |u| = 1} 

= {u G A 2 I u 2 = -1}. (76) 

Notice that we can also write every a G M x as in (74), by setting 6 = or ±7r 
according as a > or a < and u G A^ 1 ' 2 " 1 n A™ 1 * arbitrary. 

Definition 3.1. For a G A2\{0}, we define the set of periods of a by 

Period(a) := {2miru \ m G Z}, 

where u G A 2 1,2) n A% nit is the same as in (74). Recall that the pair ±u is unique 
for non-real a G A 2 \R, while u G A 2 1,2 ' n A 2 ,nit is arbitrary for a G R x . 

Notation 3.2. For a, b G A 2 , we write a = b mod(period) if b — a G Period(a). 

3.5. A multi-valued logarithmic function of an A 2 -variable. We define a 
multi-valued logarithmic function log of an A 2 -variable as the set-valued inverse 
function of the exponential function exp : A 2 — > A 2 . Explicitly, we have 

(a) if a > 0, then 

logo = {log R a + 2miru | m G Z, u G A 2 1,2) n A 2 ,nit }, (77) 
where log R : R + — > K is the usual real-valued logarithmic function; 



DELAMBRE-GAUSS FORMULAS IN HYPERBOLIC 4-SPACE 



15 



(b) if a < 0, then 

log a = {log R \a\ + (2m + l)nu \ m £ Z, u £ A^' 2) n A^ nit }; (78) 

(c) if a € A 2 \M, then 

log a = {log R |o| + (6 + 2rrm)u \ m £ Z}, (79) 

where and it are the same as in the polar decomposition (74) of a (recall that 6u 
is well dehned, as already observed in §3.4). 

3.6. A single valued logarithmic function. We may also define a single valued 
function Log : A 2 \R< — > A 2 as follows: 

(a) if a > 0, then 

Log a = log R a £R; (80) 

(c) if a € A 2 \K, then 

Log a = log M \a\ + 9u £ A 2 , (81) 

where 9 and it are the same as in the polar decomposition (74) of a; in particular, 
Ou is well defined. 

It can be easily verified that the function Log is continuous in A 2 \R<o. 

3.7. Operations and Q in A 2 . For x,y £ A 2 , wc define, as a subset of A 2 , 

x® y := log(exp(x)exp(y)). (82) 

Note that, in general, x © y ^ y © x. Observe that cosh(a; © y) and sinh(a: © y) are 
well-defined elements in A 2 and the following identities are easily verified: 

cosh(a; © y) — cosh x cosh y + sinhxsinhi/, (83) 
sinh(x © y) = sinhxeoshy + coshxsinhi/. (84) 

For convenience, we also write 

x 9 y := x © (—y) = log(exp(x)exp(— 2/)). (85) 

Since cosh(— x) = coshx* and sinh(— x) = — sinhx*, we obtain from (83) and (84): 

cosh(ir. 9 y) = cosh x cosh y* — sinh x sinh y* , (86) 
sinh(x © y) — sinh x cosh y* — cosh x sinh y* . (87) 

3.8. Mobius transformations of Aj X ' fixing both —1 and 1. The following 
proposition explains why we choose to define the hyperbolic functions cosh and sinh 
as we did in §3.3. 

Proposition 3.3. A Mobius transformation o/A 2 °' 1 ' ) fixes both —1 and 1 if and 
only if its Vahlen matrices ±A are of the form 

A — f cos ^ x sinhaA x£A (88) 

ysinha; coshay ' ~ J 

Proof. By Proposition 2.14 (with n = 2), a Mobius transformation of A^ ' 1 ^ fixes 
both —1 and 1 if and only if its Vahlen matrices ±A are of the form 

A=[ a , b ' 
b a 
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where a, b £ A2 satisfy ab* £ A 2 and aa* — bb* = 1. Since ab* = (ab*)* = 6a*, 
we have 

(a + b)(a - by = aa* - ab* + ba* - bb* = 1. 

Now let x £ log(a + b). Then a + b = exp(x). Therefore (a — b)* — exp(—a;) and 
a — b = (exp(— x))* — exp(— x*). It follows that a = coshx and b = sinhx. □ 

3.9. Mobius transformations of A^ 0,1 -* fixing both and 00. The following 
two propositions will be used in §5. They can be proved by easy calculations. 

Proposition 3.4. If a Mobius transformation o/A 2 °' 1 '' fixes each of— e2,e2,0 and 
00 then its Vahlen matrices are of the form a G IR + Me 1; \a\ = 1. 

Proposition 3.5. If a Mobius transformation o/A^ ' 1 ^ fixes each of — 1, 1, — e\, ei, 

and 00 then it is the identity transformation, that is, its Vahlen matrices are ±7. 

4. Geometry of SO (3) via the Clifford algebra A2 

In this section we study the geometry of SO(3) via the Clifford algebra A2 or 
quaternions. As explained earlier, the unit elements in A2 form the spinor group 
Spin(3), the double covering of SO (3). More precisely, every unit element a in A 2 
acts as a rotation p a on the Euclidean 3-space A^ fixing the origin by p a {x) = 
ax(a')~ 1 , and it is clear that p a = p- a . Conversely, every rotation p of A 2 °' 1 ' ) fixing 
the origin corresponds to exactly one pair of unit elements ±a in A2 such that 
p = p a = p-a- We shall first determine the axis and the rotation angle of p a in 
terms of a, and then study the Euler decomposition of a and determine the Euler 
angles associated to p a . This will be used in §5 to define the quaternion half lengths 
between two oriented line-plane flags in H 4 with a common perpendicular line. 

4.1. Axis and rotation angle of an element in SO(3). It is a well-known fact 
that every non-identity element <p £ SO(3) has 1 as a simple eigenvalue and hence 
fixes pointwise a unique straight line in A 2 0,1 \ called its axis. Thus <f> acts on the 
Euclidean space A 2 0,1 ' ) as a rotation about its axis. 

Let v £ A^ ' 1 - 1 be a unit eigenvector of 4> with eigenvalue 1 and let v 1 - denote the 
orthogonal complement plane of v in A^ ' 1 - 1 . We orient the Euclidean space A^ ' 1 ' 1 
by the ordered basis (1, e\, e-i). Then v £ A 1 ^' 1 ^ naturally induces an orientation of 
the plane v 1 - and cf>\ v ± : v 1 - — > v 1 - is a rotation of a certain angle 26 £ (0, 2ir) about 
the origin. We then say that <j> is the rotation of angle 29 about v (precisely, about 
the oriented axis with the orientation given by v). 

Observe that the rotation angle of <f> about —v is 2ir — 29. 

4.2. The unit elements in A2 as rotations of A 2 °' 1 ' ) . Recall that I^ = A2\{0} 
and SO(3) ^ r 2 /R x = Ag nlt /{±1}, where a pair ±a £ A 2 ,nit correspond to p a = 
p-a € SO(3) : A 2 0,1) ->■ A^ 04) defined by p a (x) = ax(a')" 1 = axa* . 

As a special case of the polar decomposition (74), we obtain 

Proposition 4.1. An element a £ A2 nlt \{±l} can be written uniquely in the form 
a = cos 9 + veie2 sin 9 = exp(#ueie2), (89) 
where 9 £ (0,7r) and v £ A^ 0,1 "* fl A™ 1 ' (that is, v is a unit para-vector). □ 
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Observe that, with the same 9 and v as in (89), we have 

— a = cos(7r — 9) + (— v)eie2 sin(-7r — 9). (90) 

Example 4.2. Let 9 e (0, 7r). Then we have 

(a) v = e2 if a = exp(6*ei), since e\ = e2e\e2] 

(b) v = —e± if a = exp(6*e2), since e2 = (— e.\)e.\e.i\ 

(c) t> = 1 if a = cxp(#eie2). 

The axis and rotation angle of p a are related to 9 and v in (89) as follows. 

Proposition 4.3. If a € A2 mt \{±l} is written as in (89), that is, 

a = cos9 + ve\e2 sin# = cxp(0ve-\_e2) 

wiffc 6 E (0, rr) and v € A 2 0,1) n A™' 4 , tfien u is an eigenvector of p a £ SO(3) with 
eigenvalue 1 and the rotation angle of p a about v is 28. 

Proof. Noticing that vv' = \v\ 2 — 1, we have 

p a (v) = ava* = (cos# + ve\e2 sin#)v(cos# + eie\v wiff) 

= (cos 2 9)v + (sin 2 9)ve\e2ve2e\v + (cos0sin0)i;(e2ei + eie2)v 
= (cos 2 9)v + (sin 2 6)vv' e\e2e2e\v = v. 

Thus v is an eigenvector of p a with eigenvalue 1. 

To determine the rotation angle of p a about v, let us first consider the special 
case where v = e 2 . In this case, a — a* = cos 9 + eisin# and p a restricted to 
the plane (e^ 1 - = R + Rei oriented by e2 is a rotation about the origin of angle 
29 E (0, 2tt) since we have 

p a (l) = ala* = a 2 = cos(26>) + e x sin(26>). 

In the general case, we choose c € A 2 \{0} so that p c (v) = cv(c') -1 = e 2 . Then 

cac -1 = cos^ + (sin#)cueie2C = cos^ + (sm#)cv(c') _1 eie2 = cos# + e\ sin#. 

By the special case just considered, p cac -i = p c p a (p c )~ 1 is the rotation of angle 29 
about e2- Hence p a is the rotation of angle 29 about v. □ 

Remark 4.4. It follows from (90) and Proposition 4.3 that the rotation angle of 
p_ a = p a about — v is 2(ir — 9). 

4.3. Euler decomposition of units in A2. Given a £ A™ 1 ', suppose 

a = a + a%ei + a 2 e 2 + ai 2 eie 2 , a , a%, a 2 , 012 G K. 

We pursue the following Euler decomposition: 

a = exp(aei) exp(/3eie2) exp(7ei), (91) 

with a,/3,7 € M/27rZ. (A geometric interpretation of this decomposition will be 
given in the next subsection.) Expanding the right side of (91), we obtain 

a = cos /3 cos(7 + a) + e\ cos /3 sin(7 + a) + e2 sin /? sin(7 — a) 

+ eie2sin/3cos(7 — a). (92) 
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Thus the decomposition (91) is equivalent to the following system of equations: 



a = cos /3 cos (7 + a), (93) 

ai = cos/3sin(7 + a), (94) 

a 2 = sin/3cos(7 — a), (95) 

012 = sin /3 sin(7 — a). (96) 



It follows that cos 2 /3 = a 2 , + a 2 and sin 2 /? = a\ + a\ 2 . If sin(2/3) 7^ 0, then a 
solution of /3 € R/2-7rZ determines 7 + a, 7 — a e K/27rZ and hence we obtain 
two solutions, (a,/3, 7) and (a + tt, /3, 7 + 7r), to the above system of equations in 
(R/27rZ) 3 . Since there are exactly four solutions of (3 G R/27rZ, we obtain in total 
eight solutions to the system of equations in (K/27rZ) 3 as follows: 

(a, /3, 7), (a + 7r,^,7 + 7r); 

(a + 7r,/3 + 7r,7), (a,/3 + vr,7 + 7r); 

(a + §7r, -/3, 7 - |7r), (a - ±tt, -/3, 7 + |7r); 

(a — i7T, — /3 + 7T, 7 — |7r), (a + ^TT, — P + 7T, 7 + |7r). 

It is easy to verify that the condition sin 2/3 = is equivalent to that 

a e (R + Kei) U (Re 2 + Keie 2 ) = Ai U (A x e 2 ). 

In this case, sin p = or cos P = and only a + 7 € M/27rZ or only a — 7 6 R/27rZ 
is determined. As a result, there arc infinitely many solutions in (R/2-7rZ) 3 . 

Definition 4.5. We call a triple (a, /3, 7) € (R/2ttZ) 3 or (a, /3, 7) e (M/ttZ) 3 rerjwkr 
if and only if sin 2/3 7^ 0. For convenience, we introduce the following notations: 

(M/2vrZ) 3 cg := {(a,/3, 7 ) G (R/2ttZ) 3 | sin 2/3 ^ 0}; (97) 

(R/7rZ) 3 cg := {(a,p,j) € (R/ttZ) 3 | sin 2^ ^ 0}. (98) 

Then we have proved the following propositions. 

Proposition 4.6. The map $ : (K/27rZ) 3 cg A 2 lnit \(Aie 2 U A x ) de/med by 

(a, /3, 7) 1 — > a = exp(aei) exp(/3eie2) exp(7ei) 

is an eight-fold covering map. 

Proposition 4.7. The map * : (M/7rZ) 3 cg -> (A^ nit \(Aie 2 U Ai))/{±1} de/med by 

(a, /?, 7)1 — >• ±a = exp(aei) exp(/3eie 2 ) exp(7ei) 

is a two-fold covering map. Explicitly, we have 

*- x (*(a, /?, 7)) = {(«,j8,7). (« + § 7r '-A7+ 5*0}- 

4.4. Geometric interpretation of the Euler decomposition (91). Given an 
element a = exp(aei) exp(/3eie 2 ) exp(7ei) G A™ 1 * with a, /3, 7 £ R/27rZ, the or- 
thogonal transformation p a : A^ 0,1 - 1 — > A^ ' 1 ' defined by p a (x) = ax(a')" 1 can be 
decomposed as follows: 

Pa = Poxp(Qe 1 )Poxp(/3eie 2 )Pcxp(7e 1 ) = V^Wl, (99) 

where 771 = p oxp ( Qei ), 772 = mPcxp^ei^r 1 and % = (wOPexp^eO^i) -1 - 

By Proposition 4.3, the orthogonal transformation 771 = p CX p( aei ) is the rotation 
of angle 2a about e 2 , the orthogonal transformation r? 2 = '7iPcxp(/3eie 2 ) 7 7r 1 ' s the 
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rotation of angle 2/3 about 771 (1) = cxp(2aei), and the orthogonal transformation 
r]3 = (wOPexpfreijCwi) -1 is the rotation of angle 27 about 772771(62) = 772(62). 

The transformation p a — 773772771 G SO(3) can be better understood if we consider 
its action on the unit sphere § 2 in the Euclidean space A^ 0,1 "*: 

§ 2 = A^ 04) n Aj 1 "* = {x + x x ex + x 2 e 2 \ x , X!,x 2 e M, x% + x\ + x\ = 1}. 

Note that SO(3) acts freely and transitively on the space of ordered orthogonal 
frames of A^ ' 1 ^ and hence on the unit tangent bundle of § 2 : 



Thus an element </> G SO(3) is determined by the image of a chosen unit tangent 
vector, say (l,ei), of S 2 under <fi. Alternatively, we may think of a unit tangent 
vector (x, u) of § 2 as a pointed, oriented great circle C x _ u in S 2 since the unit vector 
u G §> 2 is tangent at x G § 2 to a unique great circle C x ^ u passing through x and 
gives C x ^ u an orientation. 

Given a = exp(aei) exp(/3eie2) exp(7ei) € A™ 1 ', recall that p a = ri^rj 2 rii as in 
(99). Now we explain how to obtain geometrically the unit tangent vector 



For this, we start from (1, ei) G T 1 S 2 , traverse along the oriented great circle C\^ ei 
in its orientation by distance 2a G [0, 2ir) and arrive at 



Then we perform the rotation about 771 (ei) by angle 2/3 G [0, 2ir) and arrive at 



Finally, we traverse along the the oriented great circle C mm ^i) , n2rll ( ei ) m its orien- 
tation by distance 27 G [0, 2n) and arrive at (p a (l), p a (ei)) G T 1 ^ 2 . 



T 1 ^ 2 = {(x,u) G § 2 x § 2 I (x,u) = 0}. 



p a (l,ei) = (p„(l),p„(ei)) eT'S 2 . 



7 ?1 (l,e 1 ) = (7 ?1 (l),77 1 (e 1 ))GT 1 § 2 . 



f?277i(l,ei) = (772771(1), 772771 (ei)) G T^ 2 . 




Figure 2. The Euler angles 2a, 2/3, 27 
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Conversely, given a general unit tangent vector (a;, u) £ T 1 ^ 2 , we may read off a 
pair of elements ±et £ A™ 1 ' as follows so that the transformation p a = p- a sends 
(1, ei) to (x, u). For this, consider the pointed, oriented great circles C\^ ei and C XM . 

First suppose these two great circles are distinct. Let z be one of the intersection 
points of them (the other one is —z). Suppose that, starting from (l,ei) £ T^S 2 , 
we traverse along the oriented great circle Ci iei in its orientation by distance 2a £ 
[0, 2ir) to arrive (z, w) £ T 1 ^ 2 . Suppose further that we need to perform a rotation 
by angle 2(3 £ [0, 2ir) about z to arrive at (z, v) £ T 1 § 2 so that the great circles 
C ZtV and C XtU coincide and are orientated in the same way. Finally, let 27 £ [0, 2n) 
be the angle traversed from z to x along the oriented great circle C z , v — C-x,u in its 
orientation. In this way, we have obtained the desired decomposition (91), that is, 
a = exp(aei) exp(/3eie2) exp(7ei). See Figure 2 for an illustration of the procedure 
described above. Note that if we choose — z instead of z, then we obtain the Euler 
decomposition (91) for —a instead of a, that is, 

- a = exp((a + ^Tr)ei) exp((7r — /3)eie 2 ) exp((7 + 5 7r)ei). (100) 

Remark 4.8. The two triples of angles (2a, 2/3, 27) and (2a + ir, 2ir - 2/3, 27 + n) 
in (R/2-7rZ) 3 are (the two choices of) the Euler angles associated to the orthogonal 
transformation p a = p_ a £ SO (3). Note that Euler angles have been widely used 
by physicists and astronomers (see for example [10, pp.289] or [17, pp.100]). 

The case when the great circles Ci, ei and C X:U coincide is simpler. In this case, 
we have either a — exp(7ei) or a = exp(-|eie2) exp(7ei) for some 7 £ [0, 7r). 

4.5. An identity of Arnold on quaternions. As a special case of (92), we obtain 
identity (101) below which was used in an essential way by Arnold in [ ]. 

Proposition 4.9 (Arnold). For s,t 6l, the following identity holds: 

exp(sei) exp(teie2) exp(— sex) = exp(£ exp(2sei)eie2). (101) 
Corollary 4.10. For a, ft, 7 £ K, there hold 

exp(aei) exp(/3eie 2 ) exp(7d) = exp(2/3 exp(2aei)eie2) exp((a + 7)ei) (102) 

= exp((a + 7)ei) exp(2/3exp(— 2ae\)e\ei). (103) 

5. Half-distances between oriented lines or line-plane flags in H 4 

In this section we define the {ei, e 2 }-quaternion half distances from one oriented 
line-plane flag (an oriented line and an oriented plane containing it) in H to another 
along an oriented line orthogonal to the two flags (when such a line exists), and the 
e\- and e2-complex half distances from one oriented line in H 4 to another along an 
oriented line-plane flag orthogonal to the two lines. 

5.1. Lines and planes in H 4 . Recall that we set the upper half-space model of 
the hyperbolic 4-space as 

H 4 = {x Q + xi€i + x 2 e 2 + x 3 e 3 \ x Q , xi,x 2 , £3 £ K, x 3 > 0} C A3 1 . 

Its conformal boundary at infinity is 9H 4 = A^ ' 1 ' = A^'^Ujoo}. In this model, an 
orientation-preserving isometry of H 4 , r\ £ Isom + (H 4 ), is exactly a Mobius trans- 
formation of A^ 0,1 - 1 given by a pair of Vahlen matrices ±A £ SL(2,T2). It is a 
well-known fact that the group Isom + (H 4 ) of orientation-preserving isometries of 
H 4 acts on the set of ordered orthogonal frames of H 4 freely and transitively. 
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Convention 5.1. We shall always orient H 4 by choosing the orientation of H 4 
given by the ordered orthogonal frame at e-$ 

Frame e3 (1, ei, e2, oo) := (ij—iji], -£<[_ ei!ei ], i[_ e2)62 ], i[o,oo] 

). (104) 

We introduce notations for oriented lines and planes in H as follows. 

Notation 5.2. For distinct points u,v G <9H 4 , we use Li u>v ] to denote the unique 
line in H 4 connecting them, oriented from u to v. 

Notation 5.3 (The horizontal and vertical lines in H 4 ). We denote respectively 
by Lh and L v the horizontal and vertical oriented lines in H : 

£h = £[-i,i], L v = L[o,oo]- (105) 

Notation 5.4. For intersecting (distinct) oriented lines L\ and L 2 in H 4 , we use 
L\ V L2 to denote the oriented plane in H 4 which contains both L\ and L 2 , with 
the orientation determined by the ordered frame (L 1; L 2 ). ln particular, when 

L\ = £[u,u] and L 2 = L[ xy ] 
with distinct u,v,x,y G <9H 4 , we shall simply write L\ V Z 2 as Ilr^iyUj,]. Note 
that, with this notation, IL[ UjV ] V [ Xj y] — IL[y )X ]v[u,v] = ^-[x,y]v[v,u] — TL[v,u]v[ y ,x] ■ 

Notation 5.5 (The horizontal and vertical planes in H 4 ). We denote respectively 
by Ilh and II V the horizontal and vertical oriented planes in H 4 : 

Ilh = U[-i j i]y[- ei ,e 1 ], II V = II[_ e2ie2 ]v[0,oo]- (106) 

5.2. Orientation-preserving isometries of H 4 acting along Lr 0)OO ]. We say 
that an isometry 77 g Isom + (H 4 ) acts along the line i[o,oo] if r\ fixes both and 00 
on 9H 4 . In this case, 77 can be decomposed as the composite, in either order, of a 
pure hyperbolic translation along £[ 0jO o] an< ^ a rotation about Xro j00 i. 
Indeed, the Vahlen matrices ±A of 77 arc of the form 

A=(l , a G T 2 = A 2 X = A 2 \{0}. 



v a* 

For A > 0, let t\ € Isom + (H 4 ) be the isometry whose Vahlen matrices are 
±^ q A -1 ) ' T^ 011 T \ is a pure hyperbolic translation along the oriented line £[0,00] 
by the signed translation distance 2 log R Ael. 

By Proposition 4.3, the element p a /\ a \ & SO(3), when regarded as in Isom + (H 4 ), 
is a rotation about a plane in H 4 which contains the line Lro,<x>] ■ 

It is easy to verify that 77 = T\ a \ p a /\ a \ = p a /\ a \ T\ a \. 

5.3. Line-plane flags in H 4 . We shall consider the geometric configuration in H 4 
which is a plane with a line in it singled out. 

Definition 5.6. A flag or line-plane flag F in H 4 is an ordered pair [L, II) where 
L is a line and II a plane in H such that L is contained in II. An oriented flag 
F = (L, II) is obtained from a flag F = (L, II) by orienting each of L and II. 

Notation 5.7 (The horizontal and vertical flags in H 4 ). We denote respectively 
by Fh and F v the horizontal and vertical oriented flags in H 4 : 

F h = (L h ,II h ), F v = {L v ,ti v ). (107) 
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Definition 5.8. (i) We say that two lines or a line and a plane in H 4 are orthogonal 
to each other if they intersect and are perpendicular to each other, (ii) We say that 
a line L' and a flag (L, II) in H 4 are orthogonal to each other if L' is orthogonal to 
both L and n, that is, L' intersects L and L' is perpendicular to II. 

Definition 5.9. By a flag-line cross, (F,L'), in H 4 we mean a flag F and a line 
L' in H 4 which are orthogonal to each other. We say that a flag-line cross in H 4 is 
oriented if both the flag and the line are oriented. 

5.4. Correspondence between flag-line crosses and orthogonal frames in 

H 4 . We set up a one-to-one correspondence between the set of oriented flag-line 
crosses in H 4 and the set of ordered orthogonal frames in H 4 as follows. 

Given an ordered orthogonal frame in H 4 , that is, an ordered quadruple of ori- 
ented lines (L\, L2, £3, L4) in H 4 meeting at a point, we obtain an oriented flag-line 
cross in H 4 such that the oriented flag is (L±, L\ V L2) and the oriented line is £4. 

Conversely, given an oriented flag- line cross (F, L') in H 4 , where F = (L, II) and 
L n L' — x, we obtain an ordered orthogonal frame (L, L2, L3, L') in H 4 as follows: 
L2 is the oriented line in II passing through x and orthogonal to L such that the 
ordered pair (L, L2) determines the positive orientation of II, and Z3 is the oriented 
line in H 4 passing through x and orthogonal to all of L, L2 and L' such that the 
ordered quadruple (L, L2, L3, L') determines the positive orientation of H 4 . 

It follows from the one-to-one correspondence just established that, given two 
oriented flag-line crosses in H 4 , there exists a unique orientation-preserving isometry 
of H 4 which sends the first oriented flag-line cross to the second. 

5.5. A specific isometry of H 4 . We shall need in §5.7 the isometry 1 e Isom + (H 4 ) 
which sends the oriented flag-line cross {F y , L-^] to (i<h,-£v), or equivalently, 



In terms of orthogonal frames of H , this is equivalent to the requirement that 



The isometry t G Isom + (H 4 ) is unique and has Vahlen matrices ±K, where 



It is interesting to note that t is an involution (j? = id), with K 1 = — K. 

5.6. The quaternion half distances between oriented flags in H 4 . We define 
two quaternion half distances from one oriented flag in H 4 to another along an 
oriented common orthogonal line of the two flags (if such a line exists). Recall that 
we have identified the algebra of quaternions EI with the Clifford algebra A2. 

Definition 5.10. Suppose F\ and F2 are two oriented flags both orthogonal to an 
oriented line L in H 4 . Let l be the unique orientation-preserving isometry of H 4 such 
that i(L) = L v and l(Fi) = F&. Then there exists a unique orientation-preserving 
isometry r\ of H 4 such that t){l{L)) — l(L) and r](t,(Fi)) = 4(^2)- Since r/(L v ) = L v , 




(108) 




(109) 




(110) 
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the two quaternion half distances from F\ to F2 along L as 

(^(.Fij-Fb) = log(±a) modulo period. 

Remark 5.11. The discussions in §4.4 allow us to read off the pair ±a (or precisely, 
±a/|a|) in the above definition geometrically from the configuration (F\,L, F2). 

The following proposition recaptures Definition 5.10 and will be used later . 

Proposition 5.12. Suppose F± and F2 are two oriented flags orthogonal to an 
oriented line LinY\ 4 . Let r G Isom + (H 4 ) be such that t(L) — L and t(Fi) = F2, 
and let l £ Isom + (H 4 ) be such that l(L) = L v and l(F\) = F\ Y . Then the isometry 

ltl~ * has Vahlen matrices ±(^ &X Q^ exp( <5) ) w ^ ere ^ = ^sC^ 1 ' ^2)- O 
5.7. The e\- and e2-complex half distances between oriented lines in H 4 . 

Definition 5.13. Suppose L\ and L2 are two oriented lines orthogonal to an ori- 
ented flag F in H 4 . Let t be the unique orientation-preserving isometry of H 4 
such that l(F) = F v and i(L±) = L^. Let rj be the unique orientation-preserving 
isometry of H 4 such that t](l(F)) = l(F) and rj(b(Li)) = l(L 2 ). It follows from 

§5.2 and Proposition 3.4 that rj has Vahlen matrices ±f " ^-1) for a unique pair 



,0 a 

±a £ (R + Kei)\{0}. We define the two e\-complex half distances from L\ to L2 
along F as 

S ( S l) (L 1 ,L 2 ) =log(±a) G (R + Rei)/27reiZ. 

r 

Definition 5.14. Suppose L\ and L2 are two oriented lines both orthogonal to an 
oriented flag F in H 4 . Let t be the unique orientation-preserving isometry of H 4 
such that l(F) = F^ and l(L\) = L v . Let 77 be the unique orientation-preserving 
isometry of H 4 such that t](l(F)) = l(F) and r/(i(Li)) = £.(£2)- By Proposition 

2.14, 77 has Vahlen matrices ^ = ^ sinh 5 cosh 5 ) ^ ^ e snown m 

Proposition 5.1G below that S £ M + Ke2) for a unique pair 6,5 + ire2 ell Re2 
mod 2-7re2. We define the two e2-complex half distances from L\ to L2 along F as 

8^ %) (L X ,L\) = log(±(a + 6)) = <5,<5 + 7re 2 £ (R + Re 2 )/27re 2 Z. (Ill) 

F 

Remark 5.15. The quaternion, e\- and e2-complex half distances defined above 
are easily seen to be invariant under orientation-preserving isometries of H 4 . 

The e\- and e2-complex half distances are simply related by the R-linear map 
X ■ R + Rei — ► R + Re2 determined by x(l) = 1 and x( e i) = e 2- 

Proposition 5.16. Suppose L\ and L2 are two oriented lines orthogonal to an 
oriented flag F in H 4 . Then 

5f\L l ,L2)=x{5f\LiM))- (112) 

Proof. With no loss of generality, we may assume that F — F v and L\ = L\ x . 
Let r G Isom + (H 4 ) be such that t(F) = F and t(L x ) = L 2 . Then r has Vahlen 

matrices ±( ex ^ ? where S — S^ 1 ^ (Li , Lo)- 
\ exp(-d*)/ F ^ ' 
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Let K be the Vahlen matrix defined as in (110) and let i G Isom + (H 4 ) be the 
isometry corresponding to ±K. Then i(F) = Fj, and l{L\) = L v . 

It is easy to check that the isometry 77 := ltl^ 1 satisfies v(l(F)) = i(F) and 
rj{i,{Li)) = i{L2), with Vahlen matrices (write e = e\ + e 2 temporarily) 

( exp(5 \ _x _ , /e(cosh<5)e _1 e(sinh<5)e _1 
y exp(-<5*)/ "~ ^(sinh^e -1 e(cosh^)e _1 

= ± /coshx(<5) sinhx(<5)\ 
ysinhx(<5) coshx(^)/ 

By definition, we have Sp ) (L 1 ,L 2 ) = x(<5) = x{S [ f\Li, L 2 ))- □ 



The following proposition recaptures Definition 5.14 and will be used later . 

Proposition 5.17. Suppose L\ and Li are two oriented lines both orthogonal to an 
oriented flag F in H 4 . Let r G Isom + (H 4 ) be such that t(F) — F and t{L{) = L 2 , 
and let 1 G Isom + (H 4 ) be such that l(F) — F^ and l(L\) = L v . Then the isometry 

ltl^ 1 has Vahlen matrices ±( C ° S ^t ^J 1 *!! where S = 8^ 2 \L-\ , Lo). □ 

V cosh smn J F v ' ' 

5.8. The quaternion and the e2-complex half distances under change of 
orientations. For an oriented line L (resp. oriented plane II) in H 4 , let L~ denote 
the same line with the opposite orientation; similarly, II~, for an oriented plena n. 
For an oriented flag F = (L, II) in H 4 , let F ^ F^ and F denote the oriented 
flags (i-,3), (L,H-) and (L~ respectively. Thus F = (p-+)+-. 

The following two propositions record the changes of the e2-complex half dis- 
tances between two oriented lines and the quaternion half distances between two 
oriented flags in H 4 under a change of orientations of the oriented lines and flags 
involved and we leave it to the reader to verify the relations contained therein. 

Proposition 5.18. Suppose L\ and L2 are two oriented lines in H 4 orthogonal to 
an oriented flag F. Let 8b(LuL 2 ) = S^ 2 \Li, L2) be the C2-complex half distances 
from L\ to L2 along F. Then we have the following relations: 

5p{L2,L\) 

5^1X2) 

8p-+(Li,L 2 ) 
Sp__(L 1 ,L 2 ) 

Proposition 5.19. Suppose F\ and F 2 are two oriented flags in H 4 orthogonal to 
an oriented line L. Let 5jt(F\,F 2 ) be the quaternion half distances from F\ to F2 



= -Sp(Li,L 2 ) mod 27re 2 ; 

-* -» 7T 

= 8p(Li,L 2 ) + —e 2 mod 27re 2 ; 

= Sp(Li,L 2 ) + — e 2 mod 27re 2 ; 

= 8p{Li,L 2 ) mod 27re 2 ; 

= -5p(Li,L 2 ) mod 27re2; 

= -Sp(Li,L 2 ) mod 27re 2 . 
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along L. Then we have the following relations: 

6 z (F 2 ,Fi) = -S z (Fi,F 2 ) mod (period); 
8 z .{Fi,F 2 ) = -8 £ (F 1 ,F 2 ) mod (period); 
5 i {F 1 ,F+-) = 5 £ (Fi,F 2 )®^e 1 e 2 mod (period); 

k(Fl>F2 + ) = ^(Fi,F 2 )e|ei mod (period); 

S £ (F+- 1 F 2 ) = ^eie 2 ®5 £ (F u F 2 ) mod (period); 

%(Ff+F 2 ) = ^e 1 ®5 z (F 1 ,F 2 ) mod (period). 

6. Generalized Delambre- Gauss formulas for oriented, augmented 
right-angled hexagons in H 4 

In this section we restate and prove our generalized Delambre-Gauss formulas for 
oriented, augmented right-angled hexagons in H 4 . The proof we shall give is in 
principle the same as that for generalized Delambre-Gauss formulas for oriented 
right-angled hexagons in H 3 which we give in §7, the appendix. 

Similar to that in the case of H 3 , a right-angled hexagon in H 4 is defined to be 
a cyclic six-tuple {L n }f l=1 of lines in H 4 , with indices modulo 6, such that for all 
n = 1, • • • ,6, the adjacent lines L n and L n+ \ are orthogonal to each other. 

Definition 6.1. An augmented right-angled hexagon in H 4 is by definition a cyclic 
six-tuple {S n }n=i, with indices modulo 6, such that Si, S3, S5 are lines and S 2 , S4, Sg 
are flags in H 4 , or Si, S3, S5 are flags and S 2 , S4, Sq are lines in H 4 , and such that 
for all n = 1, • • • , 6, S n and SV1+1 are orthogonal to each other. 

Definition 6.2. An oriented, augmented right-angled hexagon {S n }^ =1 in H 4 is 
obtained from an augmented right-angled hexagon {S n }^ =1 in H by arbitrarily 
orienting each S n as S n , n= 1, • • • ,6. 

Definition 6.3. We call any choice of one of the two quaternion or e2-complex 
half distances 6g (S n —i, SVi+i) as defined in §5.6 and §5.7 a half side-length of the 
oriented, augmented right-angled hexagon {S n }n=i along S n , n = 1, 2, • • • ,6. 

Remark 6.4. It follows that if, say, Si,S 3 ,S 5 are lines and S 2 , S^Se are flags in 
H 4 , then the half side-lengths along Si, S3, S5 are {e\, e2j-quaternions (that is, 
elements in A 2 ) modulo period, and those along S2, S4, Sq are e2-complex numbers 
(that is, elements in M + M.e 2 ) modulo 2ire 2 . 

We obtain the following generalized Delambre-Gauss formulas for the quater- 
nion and the e2-complex half side-lengths of an oriented, augmented right-angled 
hexagon in H 4 which is merely a restatement of Theorem 1.1. 

Theorem 6.5. Let {«Si}^ =1 be an oriented, augmented right-angled hexagon in H 4 
with chosen quaternion and e 2 -complex half side-lengths S n along S n , n — 1, • • ■ ,6. 



2(> 



S. P. TAN, Y. L. WONG, AND Y. ZHANG 



Then the following generalized Delambre- Gauss formulas hold: 
sinh <5i cosh 62 sinh 83 + cosh 81 cosh 82 cosh S3 

= e(sinh54 cosh 65 sinh 8q + cosh £4 cosh £5 cosh Se)*] (H4) 
sinh <5i cosh 82 cosh ^3 + cosh <5i cosh 62 sinh ^3 

= e(sinh ^4 sinh 85 sinh 8% — cosh 84 sinh 65 cosh 8q)*; (H5) 
sinh 6\ sinh 62 sinh ^3 — cosh <5i sinh 62 cosh £3 

= £(sinh<5 4 cosh S 5 cosh <5 6 + cosh 84 cosh S 5 sinh 8 e )*; (H6) 
sinh 6\ sinh 82 cosh £3 — cosh 81 sinh 82 sinh £3 

= e (sinh 84 sinh (5s cosh — cosh 84 sinh (5s sinh Sq ) * , ( 1 1 7) 

wii/i e = 1 or — 1, depending on the choice of the half side-lengths {S n } n= i- 

The main idea in our proof of Theorem 6.5 is to use an identity involving six 
isometries of H 4 , namely (121), rewritten as a different identity (122) in Lemma 6.6 
below, where the entries of the Vahlen matrices of the isometries in (122) are given 
by appropriate functions of the quaternion and the e2-complex half side-lengths 
associated to the given oriented, augmented, right-angled hexagon in H 4 . 

Lemma 6.6. Given an oriented, augmented, right-angled hexagon {S n } n=1 in H 4 , 
let t„, t n , r\ n 6 Isom + (H 4 ), n = 1, • • • ,6 be determined and defined by 

T n (S n ) = S n , T n (S n -i) = SV1+1; (118) 
in ■■= (t„---t 2 t 1 )- 1 ; (119) 

Vn ■= ^n-lTV^-l: (I 20 ) 

with indices modulo 6. Then the isometries satisfy 

r 6 • • • t 2 t 1 = id; (121) 

V1V2 ■ ■ ■ m = id; (122) 
i n {S n ) = Si, L n (S n+1 ) = S 6 , n=l,3,5; (123) 
Ln(S n ) = S 6 , o n (S n+1 ) = S 1 , n = 2,4,6. (124) 

Proof. To prove (121), write t = TQT5T4T3T2T1 . Then r £ Isom + (H 4 ) and it can be 
checked that t(Sq) — Sq and t(Si) = Si by going through the following table: 

01 — > i>i — > D3 — > 03 — > D5 — > D5 — > o\, 

06 ► 02 ? <J2 ? O4 f O4 > Oq > Dq. 

Since (Se,Si) or (Si,Sq) is an oriented flag-line cross and r leaves it invariant, it 
follows that r must be the identity isometry. 

To prove (122), first note that i = t 6 = id by (121). Then 

?7l • • -776 = Tl(rf 1 T2Tl)((T2Tl) _1 T3(T2Tl)) • • • ((t5T4T3T2Ti)~ 1 T6(t 5 T4T3T2Ti)) 

= T6T5T4T3T2T1 = id. 
To verify (123) and (124), one evaluates i~ 1 (S , e) and t~ 1 (S'i) as follows: 

t-n 1 (Si) — S n , ^(Se) = S n +i, n= 1,3,5; (125) 
t ~ 1 (S ? 6 ) = 5„, t" 1 ^) = S„+i, n = 2,4,6. (126) 
This completes the proof of Lemma 6.6. □ 
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Proof of Theorem 6.5. With no loss of generality, we may assume that Si,Sa,Ss 
are lines and 52, 54, Sq are flags in H 4 . Then each S n , n — 1, 3, 5 is one of the two 
quaternion half side-lengths of {S n }n=i along 5„, and each S n , n — 2, 4, 6 is one of 
the two e2-complex half side-lengths of {S n }n=i along S n . 

By applying an orientation-preserving isometry of H 4 , we may assume that the 
oriented flag-line cross (Sq,Si) is (Fb,L v ), that is, 

5 6 = F h , 5! = L v - (127) 

Let r n E Isom + (H 4 ), n = 1, • • • ,6, with indices modulo 6, be determined by 
T n (S n ) = S n and r n (5„-i) = S n +i. Set i n = (r„ • • •r 2 ri)~ 1 and r\ n = i n -\T n L~\. 
By Lemma 6.6 we have (122) through (124). 
By (123) and (124) we have, for n = 1,3,5, 

i n-i(5„) — Si — L v , i„_i(5„_i) — Sq = F h . (128) 

Then Proposition 5.12 applies and concludes that r\ n g Isom + (H 4 ), n = 1,3,5 has 
Vahlen matrices ±A n where 

By (123) and (124) again, we have, for n = 2,4,6, 

in-l(5„) — Sq — Fh, tri-l(5„-l) — Si — L v . (130) 

Now Proposition 5.17 applies and concludes that rj n E Isom + (H 4 ), n = 2,4,6 has 
Vahlen matrices ±A n where 

A n=[^x _w J. « = 2,4,6. (131) 



\ sinh (5 n cosh <5 n 

By (122) we have 771772 ■ • • 776 = id. Hence there exists e E { — 1, 1} such that 

A 1 A 2 ---A 6 = sI, (132) 

where / is the identity 2x2 matrix, or equivalently, 

A 1 A 2 A a = eA^ 1 Aj l A^ 1 . (133) 

Working out the products of matrices on both sides of (133) and equating the 
corresponding (1,1)-, (1,2)-, (2,1)- and (2, 2) -entries, we obtain 

exp S\ cosh <5 2 exp £3 

= e(cosh<54 exp(—5g) cosh (5g + sinh ^4 exp ^5 sinh 5g)*, (134) 

— exp Si sinh ($2 exp(— S$) 

— e(sinh^ 4 exp(— 65) cosh Sq + cosh S4 exp S5 sinh Sq)* , (135) 

— exp(— SI ) sinh <5 2 exp ($3 

= e(cosh<54 exp(— <5g) sinh Sq + sinh S 4 exp <5 5 cosh Sq)* , (136) 
exp(— S{) cosh S2 exp(— S*.) 

= e(sinh S4 exp(—S*) sinh Sq + cosh^ exp<5 5 cosh Sq)* . (137) 

Now the desired formulas (114)-(117) follow from formulas (134)-(137) above by 
performing operations (134) + (137), (134) - (137), (136) + (135), and (136) - (135), 
respectively. □ 
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Remark 6.7. It is not difficult to check that the formulas (114) through (117) can 
be rewritten as the following formulas (138) through (141), respectively: 



cosh (Si e 62 © S 3 ) + cosh (Si © S 2 © S 3 ) 






= e (cosh (84 © £5 ®S 6 ) + cosh (<5 4 < 




(138) 


sinh (Si © ^2 © S3) + sinh (Si © S 2 © S 3 ) 






= e (sinh (£4 Q S* 5 Q S* 6 ) - sinh (<5 4 € 


is 5 es*))*; 


(139) 


sinh (Si Q8* 2 e8l)- sinh (<Si © S 2 © £3) 






= e (sinh (<5 4 © ^5 © <5 6 ) + sinh(c5 4 9 


)S 5 ®S 6 ))*; 


(140) 


cosh (Si ®S 2 QS^)- cosh (Si es^es^) 






= e (cosh (Si ®S 5 qSq)- cosh (<5 4 e S* 5 © <5£))*. 


(141) 



Remark 6.8. Taking the reverse involution S H> 8* , we see that the four identities 
are left invariant, with a simple shift of indices: 123456 — > 456123. 

7. Appendix: Generalized Delambre-Gauss formulas for oriented 
right-angled hexagons in H 3 

7.1. Delambre-Gauss formulas for spherical triangles. In spherical trigonom- 
etry there are for spherical triangles the important Delambre's analogies or Gauss 
formulas, (142)-(145) below, which we call the Delambre-Gauss formulas in this 
paper. These formulas were discovered by Delambre in 1807, published in 1809, 
and were subsequently discovered independently by Gauss. They play an impor- 
tant role in deriving many other important formulas in spherical trigonometry. As 
a few examples, the Napier's analogies and the law of tangents for spherical trian- 
gles follow directly, and one can derive from them the beautiful L'Huillier Theorem 
(see [11], page 44). We leave it to the reader to derive from the Delambre-Gauss 
formulas the law of cosines and the law of sines for spherical triangles. 

In the following theorem and its two corollaries we consider a spherical triangle 
in the unit sphere having side-lengths a, 6, c € (0, n) and corresponding opposite 
interior angles a, j3, 7 € (0,n). 

Theorem 7.1 (Delambre-Gauss formulas). In a spherical triangle we have 



cos |(a + b) sin ^7 = cos k(ot + (3) cos he, (142) 

sin |(a + b) sin ^7 = cos i(a — /?) sin ^c, (143) 

cos i (a — b) cos ^7 = sin i(a + (3) cos |c, (144) 

sin i(a — b) cos ^7 = sm^(a — (3) sinic. (145) 

Remark 7.2. It is easy to show that a + b > tt if and only if a + (3 > it. 
Corollary 7.3 (Napier's analogies). In a spherical triangle we have 

sin I (a — /3) / sin 5 (a + /3) = tan 5 (a— b) / tan |c, (146) 

cos \ (a — (3) I cos i (a + (3) = tan h(a + b) / tan |c, (147) 

sin |(o - 6) /sin |(a + 6) = tan | (a - /3) / cot 57, (148) 

cos±(a - b) / cos\(a + b) = tan | (a + /3) / cot -7. (149) 
Corollary 7.4 (Law of tangents). In a spherical triangle we have 

tan 5 (a - b) /tan±(a + 6) = tan |(a - /3) / tan | (a + /3) (150) 
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and two other similar formulas. 

7.2. Delambre- Gauss formulas for hyperbolic triangles and convex pla- 
nar right-angled hexagons. We may obtain the Delambre-Gauss formulas for 
hyperbolic triangles (and consequently, the Napier's analogies as well as the law of 
tangents) by simply changing all the appearance of trigonometric functions of side- 
lengths in the identities for spherical triangles over to the corresponding hyperbolic 
trigonometric functions of the same side-lengths of hyperbolic triangles. 

Theorem 7.5 (Delambre-Gauss formulas for hyperbolic triangles). For a triangle 
in the hyperbolic plane having side-lengths a,b,c > and corresponding opposite 
interior angles a,(3,j e (0,7r), the following formulas hold: 

cosh |(a + b) sin ^7 = cos j{a + f$) cosh |c, (151) 

sinh i(a + b) sin ^7 = cos \{a — (3) sinh |c, (152) 

cosh i(a — 6) cos jl — sin i(a + (3) cosh ^c, (153) 

sinh i (a — 6) cos ^7 = sini(a — /3) sinh 5 a (154) 

We may also obtain Delambre-Gauss formulas for convex right-angled hexagons 
in the hyperbolic plane as follows. 

Theorem 7.6 (Delambre-Gauss formulas for convex right-angled hexagons in the 
hyperbolic plane). For a convex right-angled hexagon in the hyperbolic plane with 
side-lengths {l n > 0}^ =l7 the following formulas hold: 

cosh + Z 3 ) sinh jh — cosh ^(Lt + 1$) cosh 5 15, (155) 

sinh i(Zi + Z 3 ) sinh |i 2 = cosh i(Z 4 — Z 6 ) sinh |/ 5 , (156) 

cosh i(Zi — Z 3 ) cosh i/ 2 = sinh i(Z 4 + Z 6 ) cosh i/ 5 , (157) 

sinh i (/1 — Z 3 ) cosh i/ 2 = sinhi(Z 4 — Z 6 ) sinhi/ 5 . (158) 

Remark 7.7. It is easy to show that l\ < l 3 if and only if I4 < Iq. 

7.3. Delambre-Gauss formulas for oriented right-angled hexagons in H 3 . 

By definition, a right-angled hexagon in H 3 is a six-tuple {L n } e n=1 of lines in H 3 
such that, for each n modulo 6, the two lines L n and L n+ \ intersect perpendicularly. 
An oriented right-angled hexagon {L n }^ =1 in H 3 is obtained from a right-angled 
hexagon {L n }^ =1 in H 3 by orienting each side-line L n as L n for n — 1, • • • ,6. 

Definition 7.8. For an oriented right-angled hexagon {L n }^ =1 in H 3 , we define its 
complex side-length a n = (L„_i, L„ + i) along L n as follows. First, there exist 
a unique orientation-preserving isometry n of H 3 so that 

r](L n ) = L[ ,cx>] and v(L n -i) = 

Since the lines rj(L n ) and rj(L n+ i) intersect perpendicularly, n(L n+ i) is of the form 
L[- Z , z ] for some z £ C\{0}. Then we define the complex side-length by 

o~n = cr £n (L„_i,L„+i) = logz e C/2niZ. (159) 

It is clear that the complex side-lengths {er„}^ =1 are invariant under orientation- 
preserving isometries of H 3 . Furthermore, dividing a n <G C/27riZ by 2, we obtain 
two half side-lengths, 5 n and 5 n + ni, in C/27riZ. 

We obtain Delambre-Gauss formulas for oriented right-angled hexagons in H 3 . 
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Theorem 7.9 (Delambre-Gauss formulas for oriented right-angled hexagons in 
H 3 ). Given an oriented right-angled hexagon {L n }^ =1 in H 3 , let S n € C/27riZ be 
one of the two complex half side-lengths o/{L„}^ =1 along L n . Then 

cosh(5i + S3) cosh(S 2 = £ cosh(<54 + Sq) cosh(S 5 , (160) 

— sinh(<5i + S3) cosh^ = e cosh(<54 — 5$) sinh^s, (161) 

— cosh(5i — 63) sinh<5 2 = e sinh(<5 4 + Sq) cosh<5 5 , (162) 
sinh(i5i — S3) sinh^ = e sinh(<54 — Sq) sinh^s, (163) 

with e = 1 or —1, depending on the choices of the six half side-lengths {S n }^ =1 . 

Our proof of Theorem 7.9 uses an identity involving six isometries of H 3 , (165) 
in Lemma 7.10 below, associated to an oriented right-angled hexagon in H 3 . 

Lemma 7.10. Given an oriented right-angled hexagon {L n }^ =1 in H 3 , let r„ € 
Isom + (H 3 ), n = 1, ••• ,6 be determined by T n (L n ) = L n and T n (L n -i) — L n +x, 
with indices modulo 6, and set rj n = (r„_i • • • T\) T n (r n -i ■ • • Ti). Then 

r 6 ■ • • n = id, (164) 
771 ■ ■ • r/ 6 = id, (165) 
and furthermore, rj n (Li) = L\ for n — 1,3,5, and r] n (Le) = Lq for n = 2,4,6. 

Proof. To prove (164), write r = T6T5T4T3T2T1. It can be checked that t(Lq) = Lq 
and t(Li) = L\ by going through the following table: 

j T\ f T2 f T3 f Ti f T 5 f T 6 f 

h\ > L\ > L3 > L3 > L$ > L§ > h\, 

j T\ f T 2 , r Ts > T Ti T T5 i T Te i f 

J^6 f J^2 > J^2 > 1^4 > J^4 > > ^6- 

Since oriented lines Lq and L\ intersect perpendicularly and r € Isom + (H 3 ) leaves 
each of them invariant, it follows that r must be the identity isometry. 
One verifies (165) by direct calculation: 

Vl-'-Ve = T 1 {T^ 1 T 2 T 1 ){{T 2 T i y 1 T3(T 2 T 1 )) ■ ■ ■ ( (T5T4T3T2T1 ) _1 T 6 (T5T4T3T2T1 ) ) 
= T 6 T 5 T 4 T 3 T 2 T 1 = id. 

The verification of the "furthermore" part of the lemma is easy. □ 

Proof of Theorem 7.9. Since the 8 n 's are invariant under orientation-preserving 
isometries of H 3 , by applying such an isometry, we may assume that 

Le = ^[-1,1] j Li = L[ Qoo y (166) 

Then the isometry r\ n £ Isom + (H 3 ), n = 1, • • • ,6 defined in Lemma 7.10 is given 
by matrices ±A n 6 SL(2,C) where 

' exp S n 



A 



. _ / cosh S n sinh S n \ _ 
An ~ Vsinh5 n coshcW' ' ,0 - 

By (165) there holds 771 • • • rje = id. Hence there exists e € { — 1, 1} such that 

A^AsAiAsAe = el, (167) 

where / is the identity 2x2 matrix, or equivalently, 

AiA 2 A 3 = eA^A^Al 1 . (168) 
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Working out the products of matrices on both sides of (168) and equating the 
corresponding (1,1)-, (1,2)-, (2,1)- and (2, 2) -entries, we obtain 

exp 5 1 cosh 82 exp S3 

= e(cosh<54 exp(— 85) cosh^ + sinh 84 exp 85 sinhJg), (169) 

— exp<5i sinh 82 exp(— 83) 

= e (sinh 84 exp (—S 5 ) cosh S e + cosh 84 exp S 5 sinh 8 6 ), (170) 

— exp(— Si) sinh 82 exp ^3 

= £ (cosh 84 exp (—8^ ) sinh 8 6 + sinh 84 exp £5 cosh <$ 6 ), (171) 
exp(— <5i) cosh (52 exp(— 83) 

= e (sinh £4 exp (—85) sinh J 6 + cosh 84 exp £5 cosh <$ 6 ). (172) 

Now the desired formulas (160)-(163) follow from formulas (169)-(172) above by 
performing operations (172) + (169), (172) - (169), (171) + (170) and (171) - (170), 
respectively. □ 

Remark 7.11. It is interesting to observe that, by suitably changing orientations 
of the side-lines, one may obtain from the single formula (160) the other three 
formulas (161)-(163). The rule is that, if we change the orientation for only one 
L n to obtain a new oriented right-angled hexagon \L' X , ■ ■ ■ , L' 6 } and choose the 
corresponding new half side-lengths 8[, ■ ■ ■ ,8' 6 as follows (with indices modulo 6): 

8 n _i = 8 n —i — ^i, S' n — —8 n , S' n+1 — <5„_|_i + 

Ki+2 = ^n+2, S' n+ 3 = S n +3, 8 n+i = 8 n+ i, (173) 

then it is easily checked that e' — e. Following this rule, one obtains formulas (161), 
(162) and (163) from (160) by simply changing the orientation of L 6 , that of L 3 , 
and those of Lq and L3, respectively. 
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